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Abstract 

The unique features (super-low background and large sensitive volume) of the 
CTF and BOREXINO set ups are used in the CAMEO project for a high sensitivity 
study of 100 Mo and 116 Cd neutrinoless 2/3 decay. Pilot measurements with 116 Cd 
and Monte Carlo simulations show that the sensitivity of the CAMEO experiment 
(in terms of the half- life limit for 0v2[3 decay) is (3-5) TO 24 yr with a 1 kg source 
of 100 Mo ( 116 Cd, 82 Se, and 150 Nd) and «10 26 yr with 65 kg of enriched 116 CdW0 4 
crystals placed in the liquid scintillator of the CTF. The last value corresponds to 
a limit on the neutrino mass of m v < 0.06 eV. Similarly with 1000 kg of 116 CdW04 
crystals located in the BOREXINO apparatus the neutrino mass limit can be pushed 
down to m v < 0.02 eV. 



1 Introduction 



Neutrinoless double beta (0i>2f3) decay is forbidden in the Standard Model (SM) since it 
violates lepton number (L) conservation. However many extensions of the SM incorporate 
L violating interactions and thus could lead to Oz/2/3 decay [p], 3- Currently, besides the 
conventional left-handed neutrino (u) exchange mechanism, there are many other possi- 
bilities to trigger this process : right-handed v exchange in left-right symmetric models; 
exchange of squarks, sneutrinos, etc. via supersymmetric (SUSY) interactions; exchange 
of leptoquarks in models with leptoquarks; exchange of excited Majorana neutrinos in 
models with composite heavy neutrinos, and so on. In that sense Oi/2/3 decay has a great 
conceptual importance due to the strong statement obtained in a gauge theory of the weak 
interaction that a non-vanishing (V2/3 decay rate requires neutrinos to be massive Ma- 
jorana particles, independently of which mechanism induces it B. Therefore, at present 
Oz/2/5 decay is considered as a powerful method to test new physical effects beyond the 
SM, while absence of this process - established at the present level of sensitivity - would 
yield strong restrictions on parameters of manifold extensions of the SM and narrow the 
wide choice of the theoretical models. At the same time 0u2f3 decay is very important in 
light of the current status of neutrino physics (see |||). Indeed, the solar neutrino problem 
(in particular lack of 7 Be neutrinos) and the deficit of the atmospheric muon neutrino flux 
|| could be explained by means of neutrino oscillations, which require nonzero neutrino 
masses. Also indication for v^j v e oscillations was found by the LSND collaboration [|. g . 
Oscillation experiments are only sensitive to neutrino mass difference, while measuring 
the Oz/2/9 decay rate can give the absolute value of the v mass scale, and hence provide a 
crucial test of neutrino mass models. 

Despite the numerous efforts to observe Oi/2/? decay beginning from 1948 up to the 
present |I[ this process still remains unobserved. The highest half-life limits were set 
in direct experiments with several nuclides: Ti^Oz/) > 10 22 yr for 82 Se 0, 100 Mo [Q; 
T 1/2 (0z/) > 10 23 yr for 116 Cd @, 0, 130 Te and 136 Xe @; and T 1/2 (0z/) > 10 25 yr for 



76 Ge|lJ,g. 

The present theoretical and experimental status of 2(3 decay investigations [[l], § makes 
it necessary to extend the number of candidate nuclides studied at a sensitivity comparable 
to or better than that for 76 Ge (neutrino mass limit m v < 0.1 — 0.5 eV). With this aim 
we consider in the present paper the use of the super-low background liquid scintillation 
detector - the BOREXINO Counting Test Facility (CTF) @ - for high sensitivity 2/3 
decay research. 



2 The CTF and choice of candidate nuclei 



The full description of the CTF and its performance have been published elsewhere [15 



16[ |17|| . Here we recall the main features of this apparatus, which are important in the 
following. 

The CTF is installed in the Gran Sasso Underground Laboratory and consists of an 
external cylindrical water tank (011x10 m; ^1000 t of water) serving as passive shielding 
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for 4.8 m 3 of liquid scintillator contained in an inner spherical vessel of 02.1 m. High 
purity water is supplied by the BOREXINO water plant, which provides its radio-purity 
level of « 10" 14 g/g (U, Th), « 1(T 10 g/g (K natural) and <5 //Bq/Zfor 222 Rn 



The liquid scintillator is a binary solution of 1.5 g/l of PPO in pseudocumene. The 
fluorescence peak emission is at 365 nm and the yield of emitted photons is ^10 4 per 



MeV of energy deposited. The attenuation length is larger than 5 m above 380 nm [19 . 
The principal scintillator decay time was measured as ~3.5 ns in a small volume, while 
as 4.5-5.0 ns with a source placed in the center of the CTF. The scintillator is purified 
by recirculating it from the inner vessel through a Rn stripping tower, a water extraction 
unit, a Si-Gel column extraction unit, and a vacuum distillation unit. It ensures that 
232 Th and 238 U contaminations in the liquid scintillator are less than (2-5) • 10~ 16 g/g. 

The inner vessel for the liquid scintillator is made of nylon film, 500 /zm thick, with 
excellent optical clarity at 350-500 nm, which allows collection of scintillation light with 
the help of 100 phototubes (PMT) fixed to a 7 m diameter support structure inside 
the water tank. The PMTs are 8" Thorn EMI 9351 tubes made of low radioactivity 
Schott 8246 glass, and characterized by high quantum efficiency (26% at 420 nm), limited 
transit time spread (cr = 1 ns), good pulse height resolution for single photoelectron pulses 
(Peak/Valley = 2.5), low dark noise rate (0.5 kHz), low after pulse probability (2.5%), 
and a gain of 10 7 . The PMTs are fitted with light concentrators 57 cm long and 50 cm 
diameter aperture. They provide 20% optical coverage. The number of photoelectrons 
per MeV measured experimentally is (300 ± 30)/MeV on average. An upgrade of the 
CTF was realized in 1999 when an additional nylon barrier against radon convection and 
a muon veto system were installed. 

For each event the charge and timing (precision of 1 ns) of hit PMTs are recorded. 
Each channel is doubled by an auxiliary channel to record all other events coming within 
a time window of 8 ms after the first event. For longer delay the computer clock is used 
(accuracy of ~0.1 s). Event parameters measured in the CTF include: 

- the total charge collected by the PMTs during 0-500 ns (event energy); 

- the tail charge (48-548 ns) used for pulse shape discrimination; 

- PMT timing to reconstruct the event in space (resolution of 10-15 cm); 

- the time elapsed between sequential events, used to tag time-correlated events. 
Due to all these measures the CTF is the best super-low background scintillator of 

large volume at present. Indeed, the total background rate in the energy region 250 - 
800 keV (so called "solar neutrino energy window") is about 0.3 count s/yr-keV- kg and 
is dominated by external background from Rn in the shielding water (^30 mBq/m 3 in 
the region surrounding the inner vessel), while the internal background is less than 0.01 
counts/yr-keV-kg. Therefore one can conclude that the CTF is the ideal apparatus for 
super-low background 2f3 decay research. 

For the choice of 2(3 candidate nuclei let us express the Oz/2/9 decay probability (ne- 
glecting right-handed contributions) as follows [§]: 



Tlj 2 ) 1 = CfiU ■ \ME\ 2 ■ (m v ) 2 
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where \ME\ is the nuclear matrix element of the Qi>2(3 decay, and G^ m (Z,Qpp) is the 
phase space factor. Ignoring for the moment the \ME\ calculation 0, it is evident from 
Eq. (1) that for the sensitivity of the 2(3 decay study with a particular candidate the 
most important parameter is the available energy release (Qpp)- First, because the phase 
space integral (hence, Oi/2/5 decay rate) strongly depends on Qpp value (roughly as Qpp). 
Second, the larger the 2(3 decay energy, the simpler - from an experimental point of view 
- to overcome background problems. The crucial value is 2614 keV which is the energy 
of the most dangerous 7's from 208 T1 decay ( 232 Th family). Among 35 candidates there 
are only six nuclei with Qpp larger than 2.6 MeV P0| : 48 Ca (Qpp = 4272 keV, natural 
abundance 5 = 0.187%), 82 Se (Qpp = 2995 keV, 5 = 8.73%), 96 Zr (Qpp = 3350 keV, 
S = 2.80%), 100 Mo (Qpp = 3034 keV, S = 9.63%), 116 Cd (Qpp = 2805 keV, 5 = 7.49%), 
and 150 Nd (Qpp = 3367 keV, 5 = 5.64%). The values of the phase space integral G^ m for 
these candidates are (in units 10~ 14 yr): 6.4 ( 48 Ca), 2.8 ( 82 Se), 5.7 ( 96 Zr), 4.6 ( 100 Mo), 4.9 
( 116 Cd), and ^20 ( 150 Nd) ||. In comparison, G^ m for 76 Ge is equal ^0.6 (in the same 
units) because of the lower 2(3 decay energy (Qpp = 2039 keV). From this list 100 Mo and 
116 Cd were chosen as candidates for 2(3 decay study with the CTF in the first phase. The 
main reason for the choice of 100 Mo - in addition to its high Qpp value - is the fact that 
the Institute for Nuclear Research (Kiev) possesses ~1 kg of 100 Mo enriched to ~99%. 
The 116 Cd was chosen because during the last decade the INR (Kiev) has developed and 



performed 2(3 decay experiments with this nuclide [pi] , |22| , |10| , which can be considered 
as a pilot step for the proposed project with the CTF. 



3 CAMEO-I experiment with luu Mo in the CTF 

There are two different classes of 2(3 decay experiments: (a) an "active" source experi- 
ment, in which the detector (containing 2(3 candidate nuclei) serves as source and detector 
simultaneously; (b) an experiment with "passive" source which is introduced in the de- 
tector system j]]. The sensitivity of any 2(3 decay apparatus is determined first, by the 
available source strengths (mass of the source), and second, by the detector background. 
Another factor very essential for determining the sensitivity is the energy resolution of 
the detector. Indeed, for a detector with poor energy resolution the events from the high 
energy tail of the 2v2(3 decay distribution run into the energy window of the Qv2(3 decay 
peak, and therefore generate background which cannot be discriminated from the 0z/2/3 
decay signal, even in principle. All of the decay features are similar: the same two parti- 
cles are emitted simultaneously from one point of the source, in the same energy region 
and with identical angular distribution. However, the better the energy resolution - the 
smaller the 2v tail becomes within the Ou interval, and thus the irreducible background 
becomes lower. Hence, we can conclude that the ultimate sensitivity to detect 0z/2/3 de- 
cay is really limited by the energy resolution of the detector, which is the most crucial 
parameter for any kind of set up for 2(3 decay study. 

For the second class of experiments the sensitivity is also restricted by the trade-off 
between source strengths and detection efficiency. The number of 2(3 decay candidate 
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nuclei can be enlarged by increasing the source thickness, which at the same time will 
lead to lower detection efficiency caused by absorption of electrons in the source and 
transformation of the measured 2(3 decay spectra (broadening of the peak and shifting it 
to lower energies). 

These experimental considerations are illustrated in Fig. 1, where results of model 
experiments to study 2(3 decay of 100 Mo are presented. The following assumptions were 
accepted for simulation: mass of 100 Mo source is 1 kg (f«6-10 24 nuclei of 100 Mo); measuring 
time is 5 years; half-life of 100 Mo two neutrino 2(3 decay Ti/ 2 (2v2(3) = 10 19 yr (e.g. see ref. 
P3fl), while for Oz/ mode T!/ 2 (0z/2/3) = 10 24 yr. The simulations were performed with the 
GEANT3.21 package |24| and event generator DECAY4 f25| , which describes the initial 



kinematics of the events in a, /3, and 2(3 decay (how many particles are emitted, their 
types, energies, directions and times of emission). The initial 2(3 decay spectra are shown 
in Fig. la and Fig. lb for different vertical scales. These spectra simulate 2(3 decay of 
100 Mo nuclei placed in an ideal detector ("active" source technique) with 100% efficiency 
for 2(3 decay events and with zero background (energy resolution and energy threshold of 
10 keV are presumed). For the next step the 100 Mo source was introduced in the same 
detector as a foil ("passive" source technique). The simulated spectra are depicted in 
Fig. lc (thickness of 100 Mo foil is 15 mg/cm 2 ) and Fig. Id (60 mg/cm 2 ). Finally the 
energy resolution of the detector (FWHM) was taken into account and results are shown 
in Fig. le (FWHM = 4% at 3 MeV) and Fig. If (FWHM = 8.8% at 3 MeV). It should 
be stressed that Fig. 1 represents the results of an ideal experiment, while in any real 
study the available results can only be worse by reason of the actual background, higher 
energy threshold and lower detection efficiency, etc. In fact, this is very strong statement 
because it allows to set the sensitivity limit for any real apparatus. For instance, it is 
evident from Fig. If that 0^2(3 decay of 100 Mo with half-life Ti/ 2 = 10 24 yr will hardly 
be observed by using the "passive" source technique with the following characteristics: 
(i) product of detection efficiency by number of 100 Mo nuclei «6-10 24 (e.g., 1 kg of 100 Mo 
at 100% efficiency, or 10 kg of 100 Mo at 10% efficiency); (ii) 100 Mo source thickness of 60 
mg/cm 2 ; (hi) FWHM = 8.8% at 3 MeV. At the same time, it is obvious from Fig. le 
that such a goal can be reached with similar apparatus with FWHM = 4% at 3 MeV, 
and with a 15 mg/cm 2 source. 

The last solution requires an increase of source area, which is usually limited by the 
available dimensions of the low background detectors used for 2(3 decay search. However, 
the unique features of the CTF - large sensitive volume and super-low background rate of 
the detector - permit an advanced 2(3 decay study of 100 Mo with the help of large square 
(f»7 m 2 ) and thin 100 Mo foil placed inside the liquid scintillator. 

The 100 Mo source in the CTF is a complex system placed in the inner vessel with 
liquid scintillator. It can be represented by three mutually perpendicular and crossing 
flat disks with diameter of 180 cm whose centers are aligned with the center of inner vessel 
of the CTF. Each disk is composed of three layers: inner 100 Mo source[| (thickness of ~15 
mg/cm 2 ) placed between two plastic scintillators of 1 mm thickness. The inner side of 

1 In fact, in the plane of the disk the 100 Mo source itself consists of four sectors with spacing between 
them of 12 cm, which helps in spatial reconstruction when events occur near the crossing of the disks. 
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' 'Active ' ' source-detector 




Figure 1: Simulated spectra of the model 2f3 decay experiment with 1 kg of 100 Mo. (la, 
lb) "Active" source technique: 100 Mo nuclei in a detector with 100% efficiency, zero 
background, and with 10 keV energy resolution, (lc, Id) "Passive" source technique: 
100 Mo source in the same detector with foil thickness 15 mg/cm 2 (lc) and 60 mg/cm 2 
(Id), (le) The same as (lc) but the energy resolution (FWHM) of the detector at 3 
MeV is 4%. (If) The same as (Id) but with FWHM = 8.8%. 



each plastic is coated with thin Al foil serving as light reflector, while the whole source 
" sandwich" must be encapsulated by thin transparent film (made of teflon or syndiotactic 
polypropylene) to avoid plastic dissolution in liquid scintillator. Plastic detectors can have 
a much longer decay constant (f.e. as Bicron plastic BC-444 with r « 260 ns and light 
output ~40% of anthracene) with respect to the liquid scintillator, thus their pulses can 
be discriminated easily from the liquid scintillator signals. The plastics tag electrons 
emitted from the 100 Mo source, reducing the background of this complex detector system 
significantly. The energy loss measured by the plastics are added to the electron energy 
deposit in the liquid scintillator to obtain an accurate value of the electron energy and to 
improve the energy resolution of the whole detector. 

3.1 Light collection, energy and spatial resolution 

The energy resolution of the scintillation detector depends mainly on the quality of scintil- 
lator itself, the fraction of light collected by PMTs, uniformity of light collection, quantum 
efficiency and noise of the PMT photocathodes, stability and noise of the electronics. The 
excellent liquid scintillator used in the CTF yields about 10 4 emitted photons per MeV of 
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energy deposited. In the present CTF design the actual optical coverage is 20% and the 
number of photoelectrons (p.e.) per MeV measured experimentally is (300 ± 30)/MeV 
on average. Thus with fourfold increase of light collection it would yield ^1200 p.e. per 
1 MeV or ^3600 p.e. for 3 MeV. 

To increase the actual optical coverage in the CTF, the PMTs can be mounted closer 
to the center of the detector. For instance, if 200 PMTs are fixed at diameter 5 m (and 
correspondingly the light concentrators' entrances at diameter 4 m), or 96 PMTs are 
fixed at diameter 3.8 m, the optical coverage is equal ~80%. We consider below the last 
configuration because it is the worst case for background contribution from the PMTsQ. 
Since the whole volume of the scintillator is divided by 100 Mo sources into 8 sectors, the 
PMTs are split into 8 groups of 12 PMTs each, so that one sector is viewed by one 
PMT group. Within the single sector (three mutually perpendicular reflector plates) 
scintillation photons would undergo less than 1.5 reflections on average before reaching 
the light concentrator aperture. The Monte Carlo simulation of the light propagation 



in such a geometry were performed with the help of GEANT3.21 program ||24j| . The 
emission spectrum and angular distribution of scintillation photons were added to GEANT 
code. The simulation finds that 3 MeV energy deposit would yield ^3700 photoelectrons 
allowing a measurement of the neutrinoless 2f3 decay peak of 100 Mo {Qpp = 3034 keV) 
with an energy resolution FWHM = 4%. 

This goal can be reached if the non- uniformity of light collection is corrected by using 
accurate spatial information about each event; hence, the spatial reconstruction ability 
of the CTF has to be enhanced also. The results of the Monte Carlo simulation prove 
such a possibility and show that spatial resolution of ^5-6 cm can be obtained with the 
upgraded CTF0. Primarily this is due to better light collection (increased by a factor of 
four). Secondly, it is owing to the spatial reconstruction method based on the comparison 
of pulse amplitudes from the different PMTs (within one group of 12 PMTs), and at the 
same time due to analysis of the time structure of each pulse (which can include direct 
and reflected light). 



3.2 Background simulation 

The model of the CAMEO-I set up with 100 Mo (described above) is used for the calcula- 
tions. We distinguish here between so called "/?" layers of the liquid scintillator, 15 cm 
thickQ on both sides of the complex 100 Mo source, and the rest of the liquid scintillator 
volume serving as an active shield for these main inner layers. In such a detector system 
the following energies are measured: 

2 Special R&D is in progress now to find optimal solution for the required 80% optical coverage in the 
CTF. 

3 The value obtained for the spatial resolution should be considered as indicative because in this 
preliminary phase of study a simplified model for light propagation in the CTF liquid scintillator has 
been used. 

4 These "/?" layers are separated from the total volume of the liquid scintillator by using the spatial 
information from the CTF. The thickness of 15 cm is chosen to guarantee the proper spatial reconstruction 
accuracy and the full absorption of the electrons emitted in the 2(3 decay of 100 Mo. 
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i) E pl and E\ l are the energy losses in the first and second plastic; 

E{ and El are the energy deposits in the first and second "/3" layer; 
iii) E ls is the energy loss in the liquid scintillator active shield. 

The energy threshold values of the detectors are set as E pl hr = 15 keV for the plastic 
and E\ s hr . = E^ hr = 10 keV for the liquid scintillator. The energy resolution is FWHM pl = 
VlO.8 ■ Ep 1 for the plastic and FWHM? = V4.8 ■ W for the liquid scintillator (FWHM, 
E pl and E@ are in keV). The latter corresponds to the value of 4% at 3 MeV. The following 
cuts are used in the simulation in order to recognize the double (3 decay events: 

i) E{ 1 or El 1 > Ef hr - 

ii) E{ 1 + El 1 > 300 keV; 

iii) if E? > Et hr , the corresponding Ef must be > E P ^ r , necessarily; 

iv) E ls < E\ s hr , i. e. there is no signal in the liquid scintillator active shield. 

The simulation of the background and decays of radioactive nuclides in the installation 
were performed with the help of GEANT3.21 plfl and the event generator DECAY4 |25|. 



3.2.1 Two neutrino 2(3 decays of 100 Mo 

The half-life of 2v2(3 decay of 100 Mo has been already measured as ^10 19 yr (e.g. ref. 
pl |8], ^J), hence the corresponding activity of a 1 kg 100 Mo source equals «13.2 mBq. 
The response functions of the CAMEO-I set up for 2u2(3 decay of 100 Mo with Ti/2 = 
10 19 yr as well as for Qv2(3 decay with Ti/ 2 = 10 24 yr (for comparison) were simulated as 
described above, and results are depicted in Fig. 2a. The calculated values of efficiency 
for the neutrinoless channel are 80% (within energy window 2.8 - 3.15 MeV), 74% (2.85 - 
3.15 MeV), and 63.5% (2.9 - 3.15 MeV). Background due to 2v2(3 decay distribution are 
19.7 counts (2.8 - 3.15 MeV), 6.1 counts (2.85 - 3.15 MeV), and 1.3 counts (2.9 - 3.15 
MeV) for 5 years measuring period. 



3.2.2 Radioactive contamination of the 100 Mo source 

The ~1 kg sample of metallic molybdenum enriched in 100 Mo to ~99% - which has to 
be applied in the present project - was already used in the quest for 2(3 decay of 100 Mo 
to the excited states of 100 Ru f26fl. In that experiment the radioactive impurities of the 



100 Mo source by 40 K and nuclides from 232 Th and 238 U chains were measured. Only 208 T1 
(measured activity is < 0.5 mBq/kg) and 214 Bi (measured activity is (12±3) mBq/kg) 
- due to their high energy release - could generate the background in the Oz/2/3 decay 
energy window for 100 Mo. The background problem associated with 100 Mo radioactive 
contamination was carefully investigated by the NEMO collaboration, which is going to 
begin 2(3 decay measurements with psIO kg of 100 Mo |27|. It was found that for the 



NEMO-3 detector the maximum acceptable internal activities of 100 Mo are 0.3 mBq/kg 
for 214 Bi and 0.02 mBq/kg for 208 T1 ||27|| . The intensive R&D were performed by the 
NEMO collaboration with an aim to show that these severe requirements can be reached 
by using presently available physical and chemical methods for 100 Mo purification [28 



On this basis, in our calculation the 100 Mo contamination criterion for 214 Bi has been 
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Figure 2: (a) The response functions of the CTF (5 kg-yr statistics) for 2(3 decay of 
100 Mo with Ti/ 2 (2^) = 10 19 yr and Ti/^Oz/) = 10 24 yr (solid histogram). Total simulated 
contributions due to 100 Mo contamination by 214 Bi and 208 T1 (dotted line), and from 
214 Bi and 208 T1 in the PMTs (dashed histogram), (b) The response functions of the CTF 
with 65 kg of 116 CdW04 crystals (5 yr measuring period) for 2v2(3 decay of 116 Cd (Ti/ 2 
= 2.7-10 19 yr), and Qv2(3 decay with Ti/ 2 = 10 25 yr (solid histogram). The simulated 
contribution from 208 T1 in the PMTs (dotted line) and from cosmogenic 110m Ag (dashed 
histogram), (c) The response functions of the BOREXINO set up with 1 t of 116 CdW04 
crystals (10 yr measuring time) for 2v2(3 decay of 116 Cd with T\/ 2 = 2.7-10 19 yr, and 0u2(3 
decay with T X j 2 = 10 26 yr (solid histogram). 



taken as 0.3 mBq/kg, and for 208 T1 as 0.1 mBq/kg, which is only 5 times better than the 
actual activity limit in our 100 Mo sample (0.5 mBq/kg)Q. 

The results of simulation, performed as described above, are as following: i) 21 Bi 
contribution to background within the energy interval 2.9 - 3.15 MeV is 6.5 counts/yr-kg; 
ii) 208 T1 contribution is equal to 0.06 counts/yr-kg. The mentioned impurities can be 
really dangerous for the experiment. However, there exists a possibility to reduce these 
background substantially by using information on the arrival time of each event for anal- 



ysis and selection of some decay chains in Th and U families [21, [2j|. With this aim, let 
us consider the 226 Ra chain containing 214 Bi: 222 Rn (T 1/2 = 3.82 d; Q a = 5.59 MeV) -> 
218 Po (3.10 m; Q a = 6.11 MeV) -> 214 Pb (26.8 m; Q p = 1.02 MeV) -> 214 Bi (19.9 m; 



5 We have accepted the less severe and more realistic criterion for 208 T1, because it was shown that 
chemical purification of Mo is very successful concerning 226 Ra chain impurities, while for 208 T1 the 
procedure is very difficult and less effective [E8| . 
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Qp = 3.27 MeV) -> 214 Po (164.3 ^s; Q a = 7.83 MeV) -> 210 Pb. The great advantage 
of the CAMEO-I experiment is the very thin 100 Mo source (~15 mg/cm 2 ), which allows 
detection of most of a and (3 particles emitted before or after 214 Bi decay, and tags the 
latter with the help of time analysis of the measured events. Indeed, our calculation gives 
the following values of the detection efficiencies: E\ = 55% for 214 Po (a particles); e<i = 
80% for 214 Pb (P); e 3 = 37% for 218 Po (a); e 4 = 32% for 222 Rn (a). The probability to 
detect at least one of these decays ( 214 Po or 214 Pb or 218 Po or 222 Rn) can be expressed as: 

e = l-(l-ei).(l-e 2 ).(l-e 3 )-(l-e 4 ). 

By substituting in this formula the calculated efficiency values, it yields e = 96.1%, 
which means that only «4% of the 214 Bi decays would not be tagged, i.e. 214 Bi contribution 
to background can be reduced by a factor of 25 (to the value of psO.26 counts/yr-kg). The 
expected total a decay rate from 238 U and 232 Th families in the entire 100 Mo source is 
~300 decays/day, however for an area of 10 x 10 cm it is only 0.4 decays/day, which allows 
use of the chains with half-life of 26.8 and 19.9 minutes for time analysis. The simulated 
background spectrum from the internal 100 Mo source contamination by the 214 Bi and 208 T1 
is presented in Fig. 2a, where the total internal background rate in the energy interval 
2.9 - 3.15 MeV is 0.3 counts/yr-kg or ^1.5 counts for 5 years measuring period. 

3.2.3 Cosmogenic activities in 100 Mo source 

To estimate the cosmogenic activity produced in the 100 Mo foil, we used the program 
COSMO j30] which calculates the production of all radionuclides with half-lives in the 
range from 25 days to 5 million years by nucleon- induced reactions in a given target. This 
program takes into account the variation of spallation, evaporation, fission and peripheral 
reaction cross sections with nucleon energy, target and product charge and mass numbers, 
as well as the energy spectrum of cosmic ray nucleons near the Earth's surface. 

For the CAMEO-I project cosmogenic activities were calculated for 100 Mo source en- 
riched in 100 Mo to 98.5% (other Mo isotopes: 98 Mo - 0.7%, 97 Mo - 0.1%, 96 Mo - 0.2%, 
95 Mo - 0.2%, 94 Mo - 0.1%, 92 Mo - 0.2%). It was assumed 5 years exposure period and 
deactivation time of about one year in the underground laboratory. The calculation shows 
that among several nuclides with Ti/ 2 > 25 d produced in 100 Mo source only two can give 
some background in the energy window of the 100 Mo neutrinoless 2f3 decay. These are 
88 Y (Q BC =3.62 MeV; T 1/2 =107 d) and 60 Co (Q /3 =2.82 MeV; T 1/2 =5.3 yr). Fortunately 
their activities are very low (~190 decays/yr for 88 Y and ~50 decays/yr for 60 Co), thus 
the estimated background in the energy region of 2.7 - 3.2 MeV is practically negligible: 
< 0.02 counts/yr-kg from 88 Y activity and < 0.005 counts/yr-kg from 60 Co. 

3.2.4 External background 

There are several origins of the external background in the CAMEO-I experiment, for 
example, neutrons and 7 quanta from natural environmental radioactivity (mainly from 
concrete walls of the Gran Sasso Underground Laboratory), contamination of PMTs by 
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40 K and nuclides from U and Th families, Rn impurities in the shielding water, cosmic 
muons (fi showers and muon induced neutrons, inelastic scattering and capture of muons), 
etc. From all of them only 7 quanta caused by PMT contamination and by Rn impuri- 
ties in the shielding water were simulated in the present work, while others were simply 
estimated as negligible on the basis of the results of ref. |31"fl , where such origins and 



contributions for the GENIUS project were investigated carefully. 

The radioactivity values of the EMI 9351 PMT accepted for the simulation are: 0.194 
Bq/PMT ( 208 T1); 1.383 Bq/PMT ( 214 Bi); and 191 Bq/PMT ( 40 K) [|17[[33|]. Also possible 



222 Rn activity ~30 mBq/m 3 in the shielding water (in the region surrounding the inner 
vessel) is taken into account. The model of the CAMEO-I detector system described 
above was used in the calculations, but with two differences: i) "/3" layers are considered 
as liquid scintillator blocks with dimensions 10x10x10 cm 3 ; ii) the threshold of the liquid 
scintillator active shield is increased to 30 keV. The simulation performed under these 
assumptions gives the following background rate in the 0z/2/3 decay energy interval (2.9 
- 3.15 MeV): i) 0.32 counts/yr-kg due to 214 Bi in PMT; ii) practically zero rates from 
208 T1 in the PMTs and 222 Rn in the shielding water. The total simulated background 
contributions due to 214 Bi and 208 T1 contamination of the PMTs is shown in Fig. 2a also. 
After 5 years it yields ^1.6 counts in the energy window 2.9 - 3.15 MeV. Summarizing 
all background sources for 5 years of measurements, one can obtain the total number of 
^4.4 counts in the energy range 2.9 - 3.15 MeV. 

3.3 Sensitivity of the CAMEO-I experiment with 100 Mo 

The sensitivity of the proposed experiment can be expressed in the term of a lower half-life 
limit for the 0^2/? decay of 100 Mo as following: 

T l/2 >ln2-N- V -t/S, (2) 

where N is the number of 100 Mo nuclei (^6-10 24 in our case); t is the measuring time (5 
years); 77 is the detection efficiency (63.5%); and S is the number of effect's events, which 
can be excluded with a given confidence level on the basis of measured data. Thus for the 
five years CAMEO-I experiment T1/2 > (13/S*)-10 24 yr. Taking into account the expected 
background of 4.4 counts, we can accept 3-5 events as the value for S (depending on the 
method of estimating S |34|. |3"5| ) which gives Ti/ 2 > (3-5) -10 24 yr and, in accordance with 



[36], the limit on the neutrino mass < m v >< 0.5 eV. On the other hand, it is evident 
from Fig. 2a that neutrinoless 2/3 decay of 100 Mo with half-life T1/2 = 10 24 yr can certainly 
be registered: the signal (13 counts) to background (4.4 counts) ratio is approximately 
3:1. 

Similar limits T 1/2 (0z/2/3) > (3-5)-10 24 yr can be obtained by the CAMEO-I set up 
with other nuclides, 82 Se (Q Pfj = 2996 keV), 96 Zr {Q pp = 3350 keV), 116 Cd {Q pp = 2804 
keV), and 150 Nd {Qpp = 3368 keV)Q. Due to its reasonable cost 116 Cd is the preferable 

6 We do not include 48 Ca (Qpp = 4272 keV) in that list because of its very low natural abundance 
(0.187%), and hence extremely high cost of a one kg 48 Ca source. 
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second candidate after 100 Mo. Note, however, that a half-life limit of ^5-10 24 yr for 150 Nd 
would lead - on the basis of the nuclear matrix elements calculation - to a restriction 
on the neutrino mass < m v >< 0.08 eV. 



4 High sensitivity 2(3 decay study of 116 Cd with the 
CTF 

The most sensitive 0z/2/3 results are obtained by using an "active" source technique [0J. 
We recall the highest limits Ty 2 — (l-2)-10 25 yr established for 76 Ge with the help of high 
purity (HP) enriched 76 Ge detectors @ 0~JJ, and bounds Tff 2 > ~10 23 yr set for 136 Xe 
with a high pressure Xe TPC [|T^], 130 Te with TeC>2 low temperature bolometers [IT], and 



116 Cd with 116 CdW0 4 scintillators 



Continuing this line, we propose to advance the experiment with 116 CdWO"4 to the 
sensitivity level of ?»10 26 yr by exploiting the advantages of the CTF. The idea is to place 
~65 kg of enriched 116 CdW04 crystals in the liquid scintillator of the CTF, which would 
be used as a light guide and anticoincidence shield for the main 116 CdW04 detectors 
(CAMEO-II project). To prove the feasibility of this task we are considering in the 
following discussion a pilot 116 Cd experiment, and then the design concept of the present 
proposal, as well as problems concerning the light collection, energy and spatial resolution, 
background sources and sensitivity estimates of the CAMEO-II project with 116 Cd. 



4.1 The pilot iib Cd study 

Here we briefly recall the main results of 116 Cd research performed during the last decade 
by the INR (Kiev)[] in the Solotvina Underground Laboratory (in a salt mine 430 m 
underground p7|), and published elsewhere JH], |TD|. The cadmium tungstate crystal 



scintillators, enriched in Cd to 83%, were grown for research [21]]. The light output 



of this scintillator is relatively large: ~40% of Nal(Tl) The refractive index is 2.3. 



The fluorescence peak emission is at 480 nm with principal decay time of ~14 //s [39 
The density of CdW04 crystal is 7.9 g/cm 3 , and the material is non-hygroscopic and 
chemically inert. 

In the first phase of the study only one 116 CdW04 crystal (15.2 cm 3 ) was placed inside 
a veto plastic scintillator and viewed by a PMT through a light-guide 51 cm long. Outer 
passive shielding consisted of high-purity copper (5 cm), lead (23 cm) and polyethylene 
(16 cm). The background rate in the energy range 2.7-2.9 MeV (Q2/3=2805 keV) was ~0.6 
counts/yr-kg-keV f2T[. With 19175 h statistics the half-life limit for 0^2/3 decay of 116 Cd 
was set as Ti/ 2 (0v) > 3.2-10 22 yr (90% C.L.) |§, while for neutrinoless 2(3 decay with 
emission of one (Ml) or two (M2) Majorons as T 1/2 (0^M1) > 1.2-10 21 yr and T 1/2 (0uM2) 
> 2.6-10 20 yr (90% C.L.) | 



7 From 1998 this experiment was carried out by the Kiev-Firenze collaboration jljj. 
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In 1998 a new set up with four 116 CdW04 crystals (total mass 339 g) was mounted in 
the Solotvina Laboratory. The enriched detectors are viewed by a special low background 
5" EMI tube (with RbCs photocathode) through one light-guide (010x55 cm), which 
is composed of two glued parts: quartz 25 cm long and plastic scintillator (Bicron BC- 
412) 30 cm long. The main detectors are surrounded by an active shield made of 15 
natural CdWC>4 crystals of large volume [^] (total mass 20.6 kg). These are viewed by 
a low background PMT (FEU- 125) through an active plastic light-guide (017x49 cm). 
The whole CdW04 array is situated within an additional active shield made of plastic 
scintillator 40x40x95 cm, thus, together with active light-guides, complete 4n active 
shielding of the main 116 CdW04 detectors is provided. The outer passive shield consists 
of high-purity copper (3-6 cm), lead (22.5-30 cm) and polyethylene (16 cm). The set up 
is isolated carefully against penetration of air which could be contaminated by radon. 

The multichannel event-by-event data acquisition is based on two IBM personal com- 
puters (PC) and a CAMAC crate. For each event the following information is stored on 
the hard disk of the first PC: the amplitude (energy), arrival time and additional tags. 
The second computer records the pulse shape of the 116 CdW04 scintillators in the energy 
range 0.25-5 MeV. It is based on a fast 12 bit ADC (AD9022) and is connected to the PC 
by a parallel digital I/O board (PC-DIO-24 from National Instruments) |3"9|| . Two PC- 
DIO-24 boards are used to link both computers and establish - with the help of proper 
software - a one-to-one correspondence between the pulse shape data recorded by the 
second computer and the information stored in the first PC. 

The energy scale and resolution of the main detector - four enriched crystals taken 
as a whole - were measured with different sources ( 22 Na, 40 K, 60 Co, 137 Cs, 207 Bi, 226 Ra, 
232 Th and 241 Am) as FWHM(keV) = -44 + V23.4 • E + 2773, where energy E is in 
keV. In particular, the energy resolution is equal to 11.5% at 1064 keV and 8.0% at 2615 
keV. Also the relative light yield and energy resolution for a particles were determined as 
following: a/(3 = 0.12 + l.l-10~ 5 E a and FWHM a (keV) = 0M3E a (E a is in keV). The 
routine calibration was carried out weekly with a 207 Bi source and once per two weeks 
with 232 Th. The dead time was monitored permanently by an LED optically connected 
to the main PMT (actual value «4.2%). 

The background spectrum measured during 4629 h in the new set up with four 
116 CdW0 4 crystals |10[ is given in Fig. 3, where old data obtained with one 116 CdW0 4 
crystal of 121 g are also shown for comparison. The background is lower in the whole en- 
ergy range, except for the (3 spectrum of 113 Cd {Qp = 316 keV) f\ In the energy region 2.5- 
3.2 MeV (location of expected Oz/2/? peak) the background rate is 0.03 counts/yr-kg-keV, 
twenty times lower than before. It was achieved first, due to new passive and active 
shielding, and secondly, as a result of the time-amplitude and pulse-shape analysis of the 
data. 

As an example of the time-amplitude technique we consider here in detail the analysis 

"1/2 

m i i - . i 2J2r>i .-ri. ..4. 

1/2 



of the following sequence of a decays from the 232 Th family: 220 Rn (Q a = 6.40 MeV, 7\ 
55.6 s) -> 216 Po (Q a = 6.91 MeV, T 1/2 = 0.145 s) -> 212 Pb. The electron equivalent energy 



i Abundance of 113 Cd in enriched 116 CdW0 4 crystals is «2% pi 
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Figure 3: Background spectrum of four enriched 116 CdW04 crystals (339 g) measured 
during 4629 h (solid histogram). The old data with one 116 CdW04 crystal (121 g; 19986 
h) normalized to 339 g and 4629 h (thin histogram). The model components: (a) 2u2j3 
decay of 116 Cd (fit value is T 1/2 (2u) = 2.6(1)-10 19 yr); (b) 40 K in the 116 CdW0 4 detectors 
(0.8 ± 0.2 mBq/kg); (c) 40 K in the shielding CdW0 4 crystals (2.1 ± 0.3 mBq/kg); (d) 
226 Ra and 232 Th in the PMTs. 



for 220 Rn a particles in 116 CdW04 is wl.2 MeV, thus all events in the energy region 0.7- 
1.8 MeV were used as triggers. Then all signals following the triggers in the time interval 
10-1000 ms (94.5% of 216 Po decays) were selected. The spectra of the 220 Rn and 216 Po a 
decays obtained in this way from data - as well as the distribution of the time intervals 
between the first and second events - are in an excellent agreement with those expected 
from a particles of 220 Rn and 216 Po ||10|| . Using these results and taking into account the 
efficiency of the time-amplitude analysis and the number of accidental coincidences (3 
pairs from 218 selected), the activity of 228 Th ( 232 Th family) inside 116 CdW04 crystals 
was determined to be 38(3) /iBq/kg JTU . 



The same technique was applied to the sequence of a decays from the 235 U family, and 
yields 5.5(14) /zBq/kg for 227 Ac impurity in the crystals. 

The pulse shape (PS) of 116 CdW0 4 events (0.25-5 MeV) is digitized by a 12-bit ADC 
and stored in 2048 channels with 50 ns/channel width. Due to different shapes of the 
scintillation signal for various kinds of sources the PS technique based on the optimal 
digital filter was developed, and clear discrimination between 7 rays and a particles was 
achieved f3~5fl . In the energy region 4.5- 6 MeV for a particles (or 0.8-1.2 MeV for 7 
quanta) numerical characteristics of the shape (shape indicator, SI) are as following: 
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S7 7 = 21.3 ± 2.0, and SI a = 32.5 ± 2.9. The PS selection technique ensures the possibil- 
ity to discriminate "illegal" events: double pulses, a events, etc., and thus suppress the 
background. For instance, PS selection of the background events, whose SI lie in the 
interval S7 7 + 2.4<r 7 < SI < SI a + 2Aa a (^90% of a events) yields the total a activity 
of 116 CdW04 crystals as 1.4(3) mBq/kg. The last value can be adjusted with activities 
determined by the time-amplitude analysis under the usual assumption that secular ra- 
dioactive equilibriums in some chains of the 232 Th and 238 U families (e.g. 230 Th — > 226 Ra 
chain) are broken. 

Since SI characterizes the full signal, it is also useful to examine the pulse front edge. 
It was found that at least 99% of "pure" 7 events (measured with calibration 232 Th source) 
satisfy the following restriction on pulse rise time : At(/is) < 1.24 - 0.5-i^ + 0.078-.E 2 , 
where E 1 is in MeV. Hence, the background pulses which do not pass this filter, were 
excluded from the residual P/j spectrum. 

The results of PS analysis of the data are presented in Fig. 4. The initial (without 
PS selection) spectrum of the 116 CdW04 scintillators in the energy region 1.2-4 MeV 
- collected during 4629 h in anticoincidence with the active shield - is depicted in Fig. 
4a, while the spectrum after PS selection of the (3/ r y events, whose SI lie in the interval 
5J 7 - 3.0a 7 < SI < S^ + 2.4ct 7 and At(/xs) < 1.24 - 0.5-£ 7 + 0.078-E 2 (containing 98% 
of (3/ j events), is shown in Fig. 4b. From these figures the background reduction due to 
pulse-shape analysis is evident. Furthermore, Fig. 4c represents the difference between 
spectra in Fig. 4a and 4b. These events, at least for energy above 2 MeV, can be produced 
by 228 Th activity from the intrinsic contamination of the 116 CdW04 crystals (measured 
by the time-amplitude analysis as described above). Indeed, two decays in the fast chain 
212 Bi (Q p = 2.25 MeV) -> 212 Po (Q a = 8.95 MeV, T 1/2 = 0.3 fjs) -> 208 Pb cannot be time 
resolved in the CdW0 4 scintillator (with an exponential decay time ~15 /xs P5| , [|3"9"||) 
and will result in one event. To determine the residual activity of 228 Th in the crystals, 
the response function of the 116 CdW0 4 detectors for the 212 Bi — > 212 Po — > 208 Pb chain 
was simulated with GEANT3.21 and event generator DECAY4. The simulated function 
is shown in Fig. 4c, from which one can see that the high energy part of the experimental 
spectrum is well reproduced (y 2 = 1.3) by the expected response for 212 Bi ^ 212 Po ^ 208 Pb 
decays^. Corresponding activity of 228 Th inside the 116 CdW04 crystals, deduced from the 
fit in the 1.9-3.7 MeV energy region, is 37(4) /iBq/kg, that is in a good agreement with 
the value determined by the time-amplitude analysis of the chain 220 Rn — > 216 Po — > 212 Pb 
38(3) /iBq/kg. 

To estimate the half-life limits for different neutrinoless 2/3 decay mode, the sim- 
ple background model was used. In fact, in the 0i>2f3 decay energy region only three 
background components (presented in fig. 3) are important: (i) external 7 background 
from U/Th contamination of the PMTs; (ii) the tail of the 2z/2/3 decay spectrum; and 
(iii) the internal background distribution expected from the 212 Bi — > 212 Po — > 208 Pb de- 
cay ( 228 Th chain). The limits for the neutrinoless mode of 2f3 decay are set as Ti/ 2 
> 0.7(2.5) -10 23 yr at 90%(68%) C.L. (for transition to the ground state of 116 Sn), while 

9 The rest of the spectrum below 1.9 MeV (Fig. 4c) can be explained as the high energy tail of the PS 
selected a particles. 
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Figure 4: Initial spectrum of 116 CdW0 4 crystals (339 g, 4629 h) without pulse-shape 
discrimination; (b) PS selected P/j events (see text); (c) the difference between spectra 
in Fig. 4a and 4b together with the fit by the response function for 212 Bi — > 212 Po — > 208 Pb 
decay chain. The fit value of 228 Th activity inside 116 CdW0 4 crystals is 37(4) /iBq/kg. 



for decays to the first 2± and second 0^ excited levels of 116 Sn as Ty 2 > 1.3(4.8)-10 22 
yr and > 0.7(2.4)-10 22 yr at 90%(68%) C.L., accordingly. For Ou decay with emission 
of one or two Majorons, the limits are: T 1/2 (0z/Ml) > 3.7(5.9)-10 21 yr and T 1/2 (0i/M2) 
> 5.9(9.4)-10 20 yr at 90%(68%) C.L. Also the half-life of 116 Cd two neutrino 2(3 decay is 
determined as T l/2 {2v2(3) = 2.6 ± 0.1(stat)±[j;I(syst)-10 19 yr |0 . 



For illustration the high energy part of the experimental spectrum of the 116 CdW04 
detectors measured during 4629 h (histogram) is shown in Fig. 5 together with the fit 
by the 2u2(3 contribution (Ti /2 =2.6-10 19 yr). The smooth curves 0u2f3Ml and 0z/2/3M2 
are 90% C.L. exclusion distributions of Ml and M2 decays of 116 Cd with T 1/2 =3.7-10 21 yr 
and Ti/2 = 5.9-10 20 yr, respectively. In the insert, the peak from 0u2(3 decay with Ti^Qv) 
= 1.0-10 22 yr and 90% C.L. exclusion (solid histogram) with Ty 2 (0v) = 7.0-10 22 yr are 
depicted for comparison. 

The following restrictions on the neutrino mass (using calculations ||36[| ) and neutrino- 
Majoron coupling constant (on the basis of calculation [[EJ) are derived from the experi- 
mental results obtained: m v < 2.6(1.4) eV and g M < 12(9.5)-10' 5 at 90%(68%) C.L. |TD|. 
It is expected that after «5 years of measurements a neutrino mass limit of m v < 1.2 
eV would be found. However further advance of this limit to the sub-eV neutrino mass 
domain is impossible without substantial sensitivity enhancement, which can be reached 
with a greater number of 116 CdW04 detectors (~65 kg) placed in the liquid scintillator 
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Figure 5: Part of the 116 CdW04 spectrum measured for 4629 h (histogram) and the fit 
by 2v2f3 contribution (Ty 2 = 2.6-10 19 yr). The smooth curves are excluded (90% C.L.) 
distributions of OvMl and 0^M2 2/5 decay of 116 Cd with T 1/2 = 3.7-10 21 yr and T 1/2 = 
5.9- 10 20 yr, correspondingly. In the insert the peak from Oz/2/9 decay with Ty 2 (0u) = 
1.0-10 22 yr is shown together with the excluded (90% C.L.) distribution (solid histogram) 
with T 1/2 (0u) = 7.0-10 22 yr. 



of the CTF. 



4.2 Design concept of the CAMEO-II project with nb Cd 

In the preliminary design concept of the CAMEO-II experiment, 24 enriched 116 CdW04 
crystals of large volume (~350 cm 3 ) are located in the liquid scintillator of the CTF and 
fixed at 0.4 m distance from the CTF center, thus homogeneously spread out on a sphere 
with diameter 0.8 m. With a mass of 2.7 kg per crystal (07x9 cm) the total 116 Cd mass 
is 20 kg (f^lO 26 nuclei of 116 Cd). It is proposed that 200 PMTs with light concentrators 
are fixed at diameter 5 m providing the total optical coverage of 80% (see footnote 3). 

The light output of CdW0 4 scintillator is (35-40)% of Nal(Tl) which yields ^1.5-10 4 
emitted photons per MeV of energy deposited. With a total light collection of ~80% 
and PMT quantum efficiency of ~25% energy resolution of several % at 1 MeV can be 
obtained. 

To justify this value a GEANT Monte Carlo simulation of the light propagation in this 
geometry was performed, which gave ^4000 photoelectrons for 2.8 MeV energy deposit. 
Therefore, with total optical coverage 80% the neutrinoless 2f3 decay peak of 116 Cd (Qpp 
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= 2805 keV) can be measured with energy resolution FWHM = 4%. The principal 
feasibility to obtain such an energy resolution with CdW0 4 crystal situated in a liquid 
has been successfully demonstrated by measurements with the help of a simple device. A 
cylindrical CdWC>4 crystal (40 mm in diameter and 30 mm in height) was fixed in the 
centre of a teflon container with inner diameter 70 mm. This was coupled on opposite sides 
with two PMTs Philips XP2412, so that the distance from each flat surface of crystal to 
the corresponding PMT's photocathode is 30 mm, while the gap between the side surface 
of the crystal and the inner surface of the container is 15 mm. The container was filled 
up with the pure and transparent paraffin oil (refractive index wl.5). Two PMTs work in 
coincidence and results of measurements with 207 Bi source are depicted in Fig. 6, where 
the spectrum obtained with the standard detector arrangement (CdW04 crystal wrapped 
by teflon diffuser and directly coupled to the PMT's photocathode with optical glue) is 
also shown for comparison. As evident from Fig. 6, a substantial (~42%) increase in 
light collection was obtained with CdWC>4 in the liquid. It resulted in improvement of 
the detector energy resolution in the whole energy region 50-3000 keV (see Fig. 7 where 
the spectra measured with 137 Cs, 60 Co and 232 Th sources are presented). It should be 
stressed that FWHM values (7.4% at 662 keV; 5.8% at 1064 keV; 5.4% at 1173 keV and 
4.3% at 2615 keV) are similar to those for Nal(Tl) crystals and have never been reached 
before with CdW0 4 crystal scintillators 
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Figure 6: The energy spectra of a 207 Bi source measured by a CdW0 4 crystal (040 x 30 
mm) for two arrangements: (a) standard, where the CdW04 crystal wrapped by teflon 
diffuser is directly coupled to the PMT's photocathode with optical glue; (b) the CdW0 4 
crystal is located in liquid and viewed by two distant PMTs (see text). 
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Figure 7: The energy spectra of 137 Cs, 60 Co and 232 Th sources measured by a CdWC>4 
crystal (040 x 30 mm) placed in liquid and viewed by two distant PMTs. 



Moreover, a strong dependence of the light collected by each individual PMT versus 
position of the emitting source in the crystal was found. Such a dependence can be 
explained by the large difference of the refraction indexes of CdWC>4 crystal (n = 2.3) 
and liquid scintillator (n' = 1.51), which leads to light redistribution between reflection 
and refraction processes due to changes of the source's position. 

General formulae for the angular distribution of the light emitted in the crystal and 
propagating in the liquid scintillator are quite cumbersome, and below we give expressions 
for some simple cases, when a CdW04 crystal with equal diameter and height (d = h = 2a) 
is placed in the center of the CTF detector and a light source is positioned on the crystal 
axis. Assuming a ratio of refraction indexes n/n' = v2 (which is close to the real value) 
and neglecting light absorption, the equation for the particular case with light source in 
the center of crystal is of the form: 

dW 1 r Icosfll 



dcosO 2V2 1 VI + cos 2 9 + ^ 

where 9 is the angle between the axis of the crystal (z-axis in our coordinate system) 
and direction of the photon in the liquid scintillator. This function is depicted in Fig. 
8a together with the distribution for another case with the light source on the bottom of 
crystal (on the axis again): 



dW 1 r I cos 1 



dcos# 2^/2 Vl + cos 2 



2®{9)l (4) 
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where 0(9) is the unit step function depending on angle 9, as following: 

0(9) - J < ^ < tt/2 

0(5) " I 0, n/2<9<n- 




Figure 8: (a) Angular distribution of the light emitted in the CdWC>4 crystal and prop- 
agating in the liquid scintillator, when a light source is on the crystal axis in the center 
and on the bottom of the crystal (solid line), (b) The same as (a) but with a light source 
positioned on an arbitrary point of the axis with coordinate z. (c) The difference between 
angular distributions for two points on the axis with coordinates z and (z + A z). arbitrary 
point of the axis with coordinate z. (c) The difference between angular distributions for 
two points on the axis with coordinates z and (z + Az). 



The angular distribution for a more general case, when the light source is posi- 
tioned on an arbitrary point of the axis with coordinate z (—a < z < a), is shown 
in Fig. 8b. When the location of the source is shifted from z to z + Az, the val- 
ues of angles 9 Z and it — 9 Z , for which dW/dcos9 changes sharply, are also changed: 
cos9 z = J2(a- \z\) 2 /(a 2 + (a- \z\) 2 ). The difference dW(z+Az)/d cos 9-dW(z)/d cos 9 
represented in Fig. 8c behaves as two peaks. The area of both peaks is equal to 
AS = a 2 A zj [a 2 + (a — |,2|) 2 ] 3 ^ 2 . Assuming that N is the full number of photoelectrons 
detected by all PMTs, and AiV is the difference of the number of photoelectrons due to 
the source shift Az, we find AiV = N-a 2 Az/ [a 2 + (a — |z|) 2 ] 3 ^ 2 . From the last expression 
one gets the formula for the spatial resolution in the CdW0 4 crystal by supposing that 
difference AN, which can be registered with 95% C.L., is equal AiV = 2\CN: 
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Substituting in Eq. (5) the crystal's dimensions a = 4 cm, and the total number of 
photoelectrons N = 4000 yields spatial resolution of ~ 4 mm in the center of the crystal 
(z = 0) and ~ 1.5 mm near the top or bottom of the crystal (z = ±a). Our GEANT 
Monte Carlo simulation proves these values and shows that, with a cylindrical CdWO"4 
crystal (07 x 9 cm) viewed by 200 PMTs, spatial resolution of 1-5 mm can be reached 
depending on the event's location and the energy deposited in the crystal (see, however, 
footnote 4). The simulated distributions of the photoelectron number among PMTs due 
to scintillations in CdW04 crystal are depicted in Fig. 9. The distance between the 
source's positions on the crystal axis is equal to Az = 2.5 mm (Fig. 9a), and Az = 1.5 
cm (Fig. 9b). 
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Figure 9: Simulated distributions of the photoelectron numbers among PMTs due to 
scintillations in CdW04 crystal (07 x 9 cm) placed in the CTF. The distance between 
source positions on the crystal axis equals Az = 2.5 mm (a), and Az = 1.5 cm (b). 



These interesting features of light collection from 116 CdWO"4 in the CTF would allow a 
reduction in the contribution from high energy 7 quanta (e.g. 208 T1) to background in the 
energy region of interest. Besides, the non- uniformity of light collection can be accurately 
corrected by using spatial information about each event in the CdW04 crystal, hence, 
helping to reach the required energy resolution of the detectors. 
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4.3 Background simulation 

The background simulations for the CAMEO-II experiment with 116 Cd were performed 
with the GEANT code and event generator DECAY4, and by using the model described 
above. 

4.3.1 Two neutrino 2(3 decays of 116 Cd 

The half-life of two neutrino 2/3 decay of 116 Cd has been measured as f»2.7-10 19 yr [41, 21 



i~0| . The response functions of the CAMEO-II set up for 2v2(3 decay of 116 Cd with Ti/ 2 
= 2.7- 10 19 yr, as well as for 0^2/3 decay with Ti/ 2 = 10 25 yr were simulated, and results 
are depicted in Fig. 2b. The calculated values of efficiency for the neutrino less channel 
are 86.1% (for energy window 2.7 - 2.9 MeV), and 75.3% (2.75 - 2.9 MeV). Background 
in the corresponding energy interval from 2u2f] decay distribution is 2.3 counts/yr (2.7 - 
2.9 MeV) or 0.29 counts/yr (2.75 - 2.9 MeV). 

4.3.2 Radioactive contamination of 116 CdW04 crystals 

The very low levels of radioactive impurities by 40 K and nuclides from natural radioactive 
chains of 232 Th and 238 U in the enriched and natural CdWO"4 crystals were demonstrated 
by the INR (Kiev) experiment [08 ]. On this basis the contamination criterion for 214 Bi 



and 208 T1 has been accepted in our calculation as ~10 /iBq/kg, which is equal to the 
actual activity value or limit determined for different samples of CdW04 crystals [pSj. 

The calculated background contribution from the sum of 208 T1 and 214 Bi activities 
is ^2000 counts/yr in the energy interval 2.7 - 2.9 MeV. However, applying the time- 
amplitude analysis with spatial resolution and pulse-shape discrimination technique de- 
veloped |21], |K]] this background rate can be reduced to ?«0.2 counts/yr or psI.O counts 
for 5 years measuring period. 

4.3.3 Cosmogenic activities in 116 CdW04 

For the CAMEO-II project cosmogenic activities in the 116 CdW04 detectors were calcu- 



lated by the program COSMO ||30|| . A 1 month exposure period on the Earth's surface 
was assumed and a deactivation time of about three years in the underground laboratory. 
Only two nuclides produce background in the energy window of the 116 Cd neutrinoless 2 ft 
decay. These are 110m Ag (Q p =3.0 MeV; T 1/2 =250 d) and 106 Ru (Qp « 40 keV; T 1/2 =374 
d) -> 106 Rh (Q p =3.5 MeV; T 1/2 =30 s). Fortunately 106 Ru activity is low and the time- 
amplitude analysis can be applied (Ti/ 2 =30 s), its estimated background is practically 
negligible: ?»0.1 counts/yr in the energy region 2.7 - 2.9 MeV, and 0.05 counts/yr (2.75 
- 2.9 MeV). The background from 110m Ag is quite large: ~23 (or ^20) counts/yr for 
the energy interval 2.7 - 2.9 MeV (2.75 - 2.9 MeV). However its contribution can be re- 
duced significantly by using spatial information because 110m Ag decays are accompanied 
by cascades of 7 quanta with energies > 600 keV, which would be absorbed in spatially 
separated parts of the detector giving an anticoincidence signature. Simulation under 
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the assumption that the 116 CdW04 crystal consists of small independent detectors with 
h — d = 1.2 cm, yields the residual background rates ~0.3 (or 0.2) counts/yr in the 
corresponding energy region 2.7 - 2.9 MeV (2.75 - 2.9 MeV). The simulated spectrum 
from the cosmogenic activity of 110m Ag is depicted in Fig. 2b. 

4.3.4 External background 

As in the previous case with 100 Mo from the various sources of external background only 
7 quanta due to PMT contamination and from Rn impurities in the shielding water were 
simulated, while others were simply estimated as negligible on the basis of the results of 
ref. [[H]. The radioactivity values of the EMI 9351 PMT accepted for the simulation are: 
0.194 Bq/PMT for 208 T1; 1.383 Bq/PMT for 214 Bi; and 191 Bq/PMT for 40 K 0, [33 . 
Also possible 222 Rn activity in the shielding water (in the region surrounding the inner 
vessel) at the level of ~30 mBq/m 3 was taken into account. The simulation performed 
under these assumptions finds that the only important contribution to the background 
in the vicinity of the 0i/2/3 decay energy is 208 T1 activity from the PMTs. The calculated 
values are ~0.8 and 0.05 counts/yr in the energy interval 2.7 - 2.9 MeV (2.75 - 2.9 MeV). 
However, with the help of spatial information available for each event occurring inside 
the 116 CdW04 crystal, these contributions can be reduced further to the level of ~0.08 
(or 0.005) counts/yr in the energy region 2.7 - 2.9 MeV (2.75 - 2.9 MeV). The simulated 
background contribution from 208 T1 contamination of the PMTs is shown in Fig. 2b. 
Summarizing all background sources gives ~3 counts/yr (0.6 counts/yr) in the energy 
interval 2.7 - 2.9 MeV (2.75 - 2.9 MeV). 

4.4 Sensitivity of the 116 Cd experiment 

As earlier we will estimate the sensitivity of the CAMEO-II experiment with the help of 
Eq. (2). Taking into account the number of 116 Cd nuclei (~10 26 ), measuring time of 5-8 
years, detection efficiency of 75%, and with expected background of 3-4 counts, one can 
obtain a half-life limit Tyz{Qu20) > 10 26 yr. On the other hand, it is evident from Fig. 2b 
that neutrinoless 2f3 decay of 116 Cd with half-life of ^10 25 yr would be clearly registered. 

It should be stressed that such a level of sensitivity for 0v2[3 decay cannot be reached 
in the presently running 2f3 decay experiments (perhaps only with 76 Ge), as well as for 



approved projects, like NEMO-3 |27| an d CUORICINO [ji^l , which are under construc- 
tion now. It was shown above that the sensitivity of the NEMO-3 set up is limited at 
the level of ^4-10 24 yr by the detection efficiency and energy resolution (see Fig. If). 
The CUORICINO project is designed to study 2(3 decay of 130 Te with the help of 60 
low temperature bolometers made of Te02 crystals (750 g mass each). Another aim of 
CUORICINO is to be a pilot step for a future CUORE project (not approved yet), which 
would consists of one thousand Te02 bolometers with total mass of 750 kg f4"2"fl . Despite 
the excellent energy resolution of these detectors (~5 keV at 2.5 MeV) the main disad- 
vantage of the cryogenic technique is its complexity, which requires the use of a lot of 
different construction materials in the apparatus. This fact, together with the lower 2f3 
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decay energy of 130 Te ((5^=2528 keV), makes it quite difficult to reach the same super- 
low level of background as obtained in experiments with semiconductor and scintillation 
detectors. For example, one can compare the background rate of 0.6 counts/yr-kg-keV at 
2.5 MeV from the current Milano experiment with twenty TeC>2 bolometers ||lT| , f42f with 
the value of 0.03 counts/yr-kg-keV in the energy region 2.5-3.2 MeV, which was reached 
in the Kiev-Firenze experiment with 116 CdW04 crystal scintillators flO ]. In that sense 
the CAMEO-II experiment has a great fundamental advantage because signaling from the 
116 CdW04 crystals to the PMTs (placed far away from crystals) is provided by light prop- 
agating in the super-low background medium of liquid scintillator, whereas cryogenic or 
semiconductor detectors must be connected with receiving modules by cables. These addi- 
tional materials (wires, insulators, etc.), whose radioactive contamination are much larger 
in comparison with TeC>2 crystals, Ge detectors or liquid scintillators, must be introduced 
in the neighborhood of the main detectors, giving rise to additional background]^]- 

Another drawback of cryogenic detectors is their low reliability. At the same time, the 
CAMEO-II experiment with 116 CdWO"4 crystals is simple and reliable, and therefore can 
run for decades without problems and with very low maintenance cost[^]. 

Moreover, the CAMEO-II project can be advanced farther by exploiting one ton of 
116 CdW0 4 detectors (^1.5-10 27 nuclei of 116 Cd) and the BOREXINO apparatus (CAMEO- 
III). With this aim 370 enriched 116 CdW0 4 crystals (2.7 kg mass of each) would be placed 
at a diameter 3.2 m in the BOREXINO liquid scintillator. The simulated response func- 
tions of such a detector system for 2t>2(3 decay of 116 Cd with Ty 2 = 2.7-10 19 yr, as well as 
for Oz/2/3 decay with T\i 2 = 10 26 yr considering a 10-year measuring period are depicted in 
Fig. 2c. Because background in BOREXINO should be even lower than in the CTF, the 
sensitivity of CAMEO-III for neutrinoless 2f3 decay of 116 Cd is estimated as Ti/ 2 > 10 27 yr, 
while Oz/2/3 decay with half-life of ^10 26 yr can be detected. This level of sensitivity can 
be compared only with that of the GENIUS project ||32| , which is under discussion now 
and intends to operate one ton of Ge (enriched in 76 Ge) semiconductor detectors placed 
in a tank (012 x 12 m) with extremely high-purity liquid nitrogen (required demands on 
its radioactive contamination are ~ 10~ 15 g/g for 40 K, 238 U, 232 Th, and 0.05 mBq/m 3 for 
222 Rn P2[ ) serving as cooling medium and shielding for the detectors simultaneously. Let 
us estimate the sensitivity of GENIUS in the same way as that for CAMEO. One ton of 
Ge detectors with enrichment ~86% (as in the current 76 Ge experiments [|14|, |H|) would 
provide f»7-10 27 nuclei of 76 Ge, thus in the optimistic case of zero background a sensitivity 
of Ty 2 > 5-10 27 yr can be reached there. 

By the aid of Eq. (1) one can obtain an expression for the neutrino mass bound derived 

from the experimental half-life limit for 0u2/3 decay as lim {m v ) = jlimT^ • G^ m ■ \ME\ 2 } 

10 There are two origins of such a background: i) radioactive contamination of the materials introduced; 
ii) external background penetrating through the slots in the detector shielding required for the connecting 
cables. 

n It should be noted that the 116 CdWC>4 crystals produced for the CAMEO-II experiment can also 
be used as cryogenic detectors with high energy resolution (4^] . In the event of a positive effect seen by 
CAMEO-II these crystals could be measured by the CUORE apparatus; in some sense both projects are 
complementary. 
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Because of the lower 2(3 decay energy the phase-space integral G^ m for 0^2/3 decay of 
76 Ge is about ten times lower than for 116 Cd. Hence, it is evident from the last equa- 
tion that neglecting the complicated problem of nuclear matrix elements calculation, the 
CAMEO-III experiment will bring at least the same restriction on the neutrino mass as 
the GENIUS project. Indeed, on the basis of the CAMEO half-life limit Ty 2 > 10 27 yr 
and using calculations |36]. [44|] one can derive a limit on the neutrino mass of ~ 0.02 
eV at 90% C.L., which is practically equal to the value ~ 0.01 eV claimed as the main 
goal of GENIUS |32|. At the same time it is obvious that the technical tasks, whose 



solutions are required for the realization of these super-high sensitive projects (GENIUS, 
CUORE, and CAMEO) are simpler for CAMEO. In fact, the super-low background ap- 
paratus needed for the last experiment is already running (this is the CTF) or under 
construction (BOREXINO), while for the CUORE and GENIUS proposals such unique 
apparati should be designed and constructed. 



5 Conclusions 

1. The unique features of the CTF and BOREXINO (super- low background and large 
sensitive volume) are used to develop a realistic, competitive, and efficient program for 
high sensitivity 2(3 decay research (CAMEO project). This program includes three natural 
steps, and each of them would bring substantial physical results: 

CAMEO-I. With a passive 1 kg source made of 100 Mo ( 116 Cd, 82 Se, 150 Nd) and located 
in the liquid scintillator of the CTF, the sensitivity (in terms of the half-life limit for 0u2j3 
decay) is (3-5)-10 24 yr. It corresponds to a bound on the neutrino mass m u < 0.1-0.3 eV, 
which is similar to or better than those of running ( 76 Ge), and future NEMO-3 ( 100 Mo) 
and CUORICINO ( 130 Te) experiments. 

CAMEO-II. With 24 enriched 116 CdW0 4 crystal scintillators (total mass of 65 kg) 
placed as "active" detectors in the liquid scintillator of the CTF the sensitivity would be 
^10 26 yr. Such a half-life limit could be obtained only by future CUORE ( 130 Te) and 
GENIUS ( 76 Ge) projects. Pilot 116 Cd research, performed by the INR (Kiev) during the 
last decade, as well as Monte Carlo simulation show the feasibility of CAMEO-II, which 
will yield a limit on the neutrino mass m v < 0.05-0.07 eV. 

CAMEO-III. By exploiting one ton of 116 CdW0 4 detectors (370 enriched 116 CdW0 4 
crystals) introduced in the BOREXINO liquid scintillator, the half-life limit can be ad- 
vanced to the level of ^10 27 yr, corresponding to a neutrino mass bound of ~ 0.02 eV. 

2. In contrast to other projects CAMEO has three principal advantages: 
i) Practical realization of the CAMEO project is simpler due to the use of already 
existing super-low background CTF or (presently under construction) BOREXINO ap- 
paratus; 
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ii) Signaling from the 116 CdW04 crystals to PMT (placed far away) is provided by 
light propagating in the high-purity medium of liquid scintillator - this allows practically 
zero background to be reached in the energy region of the Oz/2/3 decay peak; 

iii) Extreme simplicity of the technique used for 2/3 decay study leads to high relia- 
bility and low maintenance costs for the CAMEO experiment, which therefore can run 
permanently and stably for decades. 

3. Fulfillment of the CAMEO program would be a real breakthrough in the field of 
2f3 decay investigation, and will bring outstanding results for particle physics, cosmology 
and astrophysics. Discovery of neutrinoless 2/5 decay will clearly and unambiguously 
manifest new physical effects beyond the Standard Model. In the event of a null result, the 
limits obtained by the CAMEO experiments would yield strong restrictions on parameters 
of manifold extensions of the SM (neutrino mass and models; dark matter and solar 
neutrinos; right-handed contributions to weak interactions; leptoquark masses; bounds 
for parameter space of SUSY models; neutrino-Majoron coupling constant; composite 
heavy neutrinos; Lorentz invariance, etc.), which will help to advance basic theory and 
our understanding of the origin and evolution of the Universe. 
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